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Abstract: Twisted Eguchi-Kawai model is a useful tool for studying the large- A^ gauge 
f->^ I theory. It can also provide a nonperturbative formulation of the gauge theory on noncom- 

^^ ' mutative spaces. Recently it was found that the Z^ symmetry in this model, which is 

crucial for the above applications, can break spontaneously in the intermediate coupling 
^ ■ region. In this article, we study the phase structure of this model using the Monte-Carlo 

j^ . simulation. In particular, we elaborately investigate the symmetry breaking point from the 

weak coupling side. The simulation results show that we cannot take a continuum limit 

for this model. 
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1. Introduction 

The large- A^ gauge theories provide fruitful features to both phenomenology and string the- 
ory. They are simplified in the large- A^ limit while preserving essential features of QCD [l|. 
Additionally, dimensional reductions of ten-dimensional M = 1 super Yang-Mills theory 
(matrix model) are expected to provide nonperturbative formulations of superstring theory 
P, y, 1^, H], and can also be regarded as effective actions of D-branes Q. Furthermore, their 
twisted reduced versions, which we study in this article, can provide a nonperturbative 
formulation of the gauge theories on noncommutative spaces (NCYM) p. In order to 
study the nonperturbative nature of these theories, numerical simulations using lattice reg- 
ularizations are quite efficient. (Non-lattice simulations are also applicable for the reduced 
models. See references [^, |T0| for the recent progress.) 

In the large-A^ limit there is an equivalence between the gauge theory and its zero- 



dimensional reduction, which is known as Eguchi-Kawai equivalence ||lll|. Here, we consider 
the SU{N) gauge theory (YM) on Z?-dimensional periodic lattice with the Wilson's pla- 
quette action 

Sw = -PN J] J] Tr U^{x)U,{x + /i)[/t(x + i>)C/t(x), (1.1) 



where U^{x) (/i = 1,...,L)) G SU{N) are link variables and (3 is the inverse of the bare 
't Hooft coupling. In the large- A^ limit the space-time degrees of freedom can be neglected, 
and then this theory can be equivalent to a model defined on a single hyper-cube, 

Sek = -PN Y, Tr U^U^Ulul (1.2) 

which is called the Eguchi-Kawai model (EK model). The equality was shown by observing 
that the Schwinger-Dyson equations for Wilson loops (loop equations) in both theories are 
the same. In the EK model the loop equations can naively have open Wilson line terms, 
which do not exist in the original gauge theory side due to the gauge invariance. Therefore 
we need to assume that the global Z^ symmetry 

f/^ ^ e'^-U^, (1.3) 

which eliminates the non-zero expectation value of the open Wilson lines, is not broken 
spontaneously. However, soon after the discovery of the equivalence, it was found that the 



symmetry is actually broken for L* > 2 in the weak coupling region |12|. Although the 



naive EK equivalence does not hold, modifications were proposed for this issue. They are 
quenched Eguchi-Kawai model (QEK model) [0, |l3|, 14] and twisted Eguchi-Kawal model 



(TEK model) |]15[ . Historically, most of the works previously done were based on the TEK 
model because this model is theoretically interesting and numerically more practical (and 
this model describes the NCYM as mentioned before). 

In the TEK model, twisted boundary conditions are imposed and then the Z^ sym- 
metry is ensured in the weak coupling limit. It is not obvious whether the symmetry is 
broken or not in the intermediate coupling region. There is no guarantee for not violating 
the symmetry. Numerical simulations in the 1980s, however, suggested that the Z^ sym- 
metry is not broken throughout the whole coupling region. Then we have believed that 
the TEK model actually describes the large- A^ limit of the gauge theory. 

Recently some indication about the Z^ symmetry breaking was surprisingly reported 
in several context around the TEK model |16, 17|[^, 19|. The most relevant discussion 



for the present article was done by Teper and Vairinhos in |17l^ . They showed that the 
Z^ symmetry is really broken in the intermediate coupling region by the Monte-Carlo 
simulation for the D = A TEK model with the standard twist. Our work in this article 
is along this line and we mainly concentrate on investigating locations of the symmetry 
breaking from the weak coupling side in {(5, N) plane. By the Monte-Carlo simulation we 
clarify the linear behavior of critical lattice coupling 

Hi ~ L\ (1.4) 

where [5^ represents critical lattice coupling from the weak coupling side and L is the lattice 
size we have considered. This result means that the continuum limit of the planar gauge 



^In llSl [L9| a similar model with two commutative and two noncommutative dimensions were studied 
in the context of NCYM. In this case the instability of Zjv preserving vacuum was observed even in a 
perturbative calculation. This instability arises due to UV/IR mixing. 



theory cannot be described by the TEK model from the argument of the scaling behavior 
around the weak coupUng Umit. This discussion can be also applied to the NCYM case. 

This article is organized as follows. In the next section we review the TEK model 
briefly and fix our setup. In section ^ we show the numerical results for the Z^ symmetry 
breaking of the TEK model and find the scaling behavior (1.4). In section ^ we give 
the validation for the numerical result, and also discuss whether the TEK model has a 
continuum limit or not. 



2. Twisted Eguchi-Kawai model 

2.1 Action and Wilson loop 

In this study, we treat the D = 4 case. The TEK model |1^ is a matrix model defined by 
the partition function 



with the action 



Ztek = / JJ d[//xexp(-5TBx), 



Stek = -PN Y^ Z^,Tr U^U.Ulul 



(2.1) 



(2.2) 



where C/^ and dU^ {fx = 1,2,3,4) are link variables and Haar measure. The phase factors 



Zfj,u are 



Z 



flU 



exp{2TTin^,y/N) 



n 



flU 



riu^i G ^N- 



The Wilson loop operator also contains the phase Z{C) as 

WTEKiC)^Z{C){WiC)), 



(2.3) 



(2.4) 



where W{C) is the trace of the product of link variables along a contour C and Z{C) is 
the product of Z^^^s which correspond to the plaquettes in a surface whose boundary is 
C. This model is obtained by dimensional reduction of the Wilson's lattice gauge theory 
with the twisted boundary condition. With these definitions, the loop equations in the 
TEK model take the same form as those in the ordinary lattice gauge theory if the Z|^ 



symmetry, which we discuss in section |2.3| , is not broken. 

2.2 Twist prescriptions and classical solutions 

In the weak coupling limit, the path-integral is dominated by the configuration which gives 
the minimum to the action. This configuration Ufj, =Tu satisfies the 't Hooft algebra 



r^^r. 



'^UfJ,'- U'' /i; 



(2.5) 



and is called "twist-eater" . The most popular twist might be the minimal symmetric twist 
(standard twist) 

/ L L L\ 

-LOLL 

-L -L L 
\-L -L -L J 



n 



liu 



N = L^ 



(2.6) 



This twist represents L'^ lattice. In order to construct the classical solution for this twist, 
it is convenient to use the S'L(4, Z) transformation for the coordinates on T"* j20|. Using 



the SL{4, Z) transformation we can always rewrite the n^,y in the skew-diagonal form 
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flU 
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flU 



V^n^^V 
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-L J 



(2.7) 



where V is a SL(4,'Z) transformation matrix. This form makes the construction of the 
twist-eater easy. Here we define L x L "shift" matrix Sl and "clock" matrix Cl by 



Sl 



/O 1 • 


■■ 0\ 




: 1 








•• 


, Cl = 





•• 1 




\i 0... . 


•• Oy^ 





/I 



^2TTi/L 







\ 



„27ri-2/L 







,2-ni{L-l)/L 



(2.8) 



which satisfy the little 't Hooft algebra 

ClSl 



-2-Ki/L 



SlCl- 



Using these matrices, the twist-eater configuration for the skew-diagonal form ( |2.7] ) is easily 
constructed as 



Ti = Cl® 1l, 



T2 = Sl® 1l, 

r4 = 1l ® ^L- 



(2.10) 



From (|2.7| ) we can also construct the twist-eater configuration for the minimal symmetric 
twist (O) as 



ri= Cl(^1l, T2 = SlCl(^Cl, 
T^ = SlCl®Sl, r4= Sl®1l. 



(2.11) 



Although these forms are different only by the coordinate transformation, they can give 
different results except the weak coupling limit as seen in next section. 
Another kind of the twist we consider in this article is 



n 



flU 



/ rriL 
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-ml 








mL 
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—mL / 



N = mL^ 



(2.12) 



with classical solution 



Ti = Cl<S)1l^ Im, T2 = Sl<S)1l^ Im, 
Tg = 1l Cl Im, r4 = 1l ® ^L » Im- 



(2.13) 



While we write the twist using the skew-diagonal form here, we can always rewrite it in 
the symmetric form by the 5'L(4, Z) transformation. We call this twist "generic twist" in 
this article, and the minimal twists ( p.6| ) and (2.7) are particular cases (m, = 1) of the 



generic twist. As is well known, the TEK model can describe the NCYM theory |^, |^. 
Expanding the matrix model around noncommutative tori background, we can obtain 
noncommutative U{m) Yang-Mills theory on fuzzy tori. (Note that this interpretation is 
possible even at finite- A*".) Because fuzzy torus can be used as a regularization of fuzzy 
M'^, it is naively possible to give a nonperturbative formulation of the NCYM on fuzzy 
M^ by taking a suitable large-A^ limit in the TEK model. (See appendix ^ for details.) 
However, we will see later it is not the case because of the Zj^ symmetry breaking. In the 
NCYM interpretation the shift and clock matrices can be regarded as matrix realization 



of a fuzzy torus. From this point of view, twist prescription (2.12) provides YM theories 
on m-coincident four-dimensional fuzzy tori. 

2.3 Z^ symmetry 

The Z^ symmetry plays a crucial role in the Eguchi-Kawai equivalence. Generally, the YM 
theory with a periodic boundary condition has a critical size. If we shrink the volume of the 
system beyond the critical size, we encounter the center symmetry breaking, which is just 
the same as the finite temperature system. In the EK model, which is a single hyper-cubic 
model, the critical size corresponds to (3c ~ 0.19 in the lattice coupling. In the region less 
than the /3c ~ the strong coupling region - the center symmetry Zj^ is maintained. On 
the other hand, in the region larger than (5c — the weak coupling region - the symmetry is 
spontaneously broken, and then the EK equivalence does not hold. 

The TEK model avoids this problem by imposing the twisted boundary condition on 
the system instead of the periodic one. In the weak coupling limit the path integral is 
dominated by the vacuum configuration, which is twist-eater configurations, as we already 
mentioned. These configurations are invariant under global Z|^ transformation 

U^ ^ e'^^-U^, e*^" G Z^, (2.14) 

which is regarded as the U{1)^ symmetry in the large- A^ limit. As a result, Wtek{C) is 
zero if C is an open contour in the weak coupling limit. 

A key point is that the solution for this problem is obvious only at the classical level. 
That is to say, there is no guarantee to maintain the Z|^ symmetry if we take into account the 
quantum fluctuation. Going away from the weak coupling limit, the conflgurations fluctuate 
around the twist-eater. The situation can be displayed in the eigenvalue distribution of 
the link variables. In the weak coupling limit the A^ eigenvalues distribute regularly and 
uniformly on the unit circle in the complex plane, and then they are Z^ symmetric. If 
we decrease /3, the eigenvalues begin to fluctuate around the location of the twist-eater. 
If the fluctuation is not too large, the Z/, symmetric distribution is maintained. However, 
large fluctuation can make the uniform distribution shrink to a point, which corresponds to 
U^ = Ijsr configuration. In the strong coupling region the distribution is randomly uniform, 
and then the symmetry is restored. 



Although there is no guarantee to maintain the Z^ symmetry in the intermediate 
coupHng region, the 1980s numerical simulations suggested that the symmetry was unbro- 
ken. And this caused us to believe that the EK equivalence in the TEK model does hold 
throughout the whole coupling region. 

2.4 Limiting procedure 

As is well known, the scaling of the YM lattice theory behaves as /3 ~ loga^^ around 
the weak coupling limit, where a is the lattice spacing, and which is obtained by one-loop 
perturbative calculation of the renormalization group equation. If we wish to construct 
the TEK model which corresponds to the YM theory by the EK equivalence, the scaling 
of the TEK model should obey that of the YM theory. In the TEK model, the lattice size 
L relates to A^. (For the twist we consider in this article, the relation is A^ = mLp'.) Then, 
the YM system with fixed physical size / = aL can be obtained by the scaling 

/3~loga"^ ~logA. (2.15) 

In order to obtain the large-A limit with infinite volume, we should increase (3 slower than 
the scaling ( 2.15| ). If it is not the case, the system shrinks to a point. 



In the case of the NCYM, the scaling near the weak coupling limit is essentially same 
as the YM theory, that is, (3 ~ loga~^. (See appendix |A[) But if we wish to make the 
TEK model corresponding to the NCYM, there is a constraint a?L = al = fixed, which 
means that we take a scheme in which the noncommutative parameter 9 is fixed. Then, 
both the continuum limit and the infinite volume limit are simultaneously taken (double 
scaling limit). Regardless of difference of the constraint, the scaling for the NCYM we 
should take is the same as that of the ordinary YM ( 2.15| ) by the nature of the logarithm 
scaling. 

3. Zj^ symmetry breaking in the TEK model 

As mentioned in the previous section, the Zj^ symmetry breaking had not been observed 
in the older numerical simulation. However, there are several recent reports which indicate 



the symmetry breaking |16, 17, 19|. In |17|, the symmetry breaking in the D = 4 SU{N) 
TEK model was studied in the case of the standard twist up to A = 144 = 12^. The 
authors of [I^ performed the Monte-Carlo simulation starting both from a randomized 
configuration ("hot start") and from the twist-eater solution ("cold start"). In both cases 
the liff symmetry begins to break at A > 100 = 10^. At A = 144 the symmetry breaking 
and restoration patterns they observed are 



^N N N N N 

74 , 'jS , n,2 , n/l , njQ ,Pc n,^ 

j^ ■* ^jv ^ N N N N 



(A = 144, standard, hot start), 

(3.1) 

(A = 144, standard, cold start). 



where /?^ and /3^ are the first breaking point for the hot start and that for cold start, 
respectively. Note that although there is recovery of the symmetry for the cold start, the 
symmetry remains broken for the hot start. 



In this section we show the results of the numerical simulation for this symmetry 
breaking phenomena. In order to argue about the possibility of the continuum and large- 
N limiting procedure for this model, we mainly focus on the first breaking point for the 
cold start /5^, which depends on N."^ 

3.1 Simulation method 

In our simulation we use the pseudo-heatbath algorithm. The algorithm is based on pT| , 
and in each sweep over-relaxation is performed five times after multiplying SU{2) matrices. 
The number of sweeps is 0(1000) for each /?. We scan the symmetry breaking on the 
resolution of A/3 = 0.005, and then we always quote the value ±0.0025 as the error due 
to the resolution. Note that the breaking points are ambiguous because the breakdown of 
the Z^ symmetry is a first-order transition. As an order parameter for detecting the Z^ 
breakdown, we measure the Polyakov lines 



P, 



>^^ 



(3.2) 



3.2 Simulation results 

Minimal symmetric twist 

First of all we treat the minimal symmetric twist ( [2.61 ). This twist is the most standard 
and is also used in the paper ||I^ . In our study we only investigate the first Zj^ symmetry 
breaking point from weak coupling limit, that is, f3^ for this twist. (For more detailed 
information about the symmetry breaking phenomena, see fl% .) The obtained results are 

in table || and plotted in figure ||. The symmetry breaking points and patterns (Zj^ — ^ Z^ 
0i 



for N = 100; Z% ^^ Z^ for N > 100) are consistent with the results in up to A^ = 144. 
In this work we explore the simulation for larger A^. From figure |l] we can find clear linear 
dependence of /3^ on N{= L?) for A^ > 169. The fitted result in linear function using 
A^ > 169 data is 

/^f- ~ O.OOllA^ + 0.21. (3.3) 

A theoretical argument for this linear behavior is discussed in section |^. 

Minimal skew-diagonal twist 

Twists can be always transformed into the skew-diagonal form by S'L(4, Z) transformation 



as we mentioned in section ^^. As it were, the minimal symmetric twist ( |2.6D is equivalent 
to the minimal skew-diagonal twist ( p.7| ) in the weak coupling limit. However, both forms 
can represent different features by taking into account the quantum fiuctuation. Actually, 
the Z^ symmetry is already broken at A^ = 25. This fact enables us to observe the A^- 
dependence of the critical points easily. Not only is the symmetry breaking point different 
from the symmetric form, so is the breaking and restoration pattern. Figure shows the 
expectation value of the plaquette (top) and the Polyakov lines (besides the top) versus 



^Strictly speaking, the symmetry preserved in the weak coupling region is not Z^ but Z|,. However, Z|, 
is sufficient for the Eguchi-Kawai equivalence so we do not dare to distinguish them in this article. 
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P^ 


iV 
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P^ 


100 


10 


0.3525 ± 0.0025 


225 


15 


0.4575 ± 0.0025 


121 


11 


0.3625 ± 0.0025 


256 


16 


0.4875 ± 0.0025 


144 


12 


0.3775 ± 0.0025 


289 


17 


0.5275 ±0.0025 


169 


13 


0.3975 ± 0.0025 


324 


18 


0.5675 ± 0.0025 


196 


14 


0.4225 ± 0.0025 





Table 1: Critical lattice coupling from the weak coupling side /3^ for the minimal symmetric twist. 
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Figure 1: Plot of (3^ versus N for the minimal symmetric twist. Fit line is equation (3.3), which 
is obtained using N > 169 data. 



P for the cold start at N 
restoration pattern: 



100. For A^ > 100 we find the Z^ symmetry breaking and 



Z^ 
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73 



72 






Z^ (A^ = 100, minimal skew-diagonal, cold start), (3.4) 



which represents a difference from the symmetric form case (p.l[). The first breaking 



pattern Zj^ — ^ 



^j(f is, however, the same as that in the symmetric twist. (We note 



that for A^ < 81 the first breaking pattern is Z^ 
" '^■^ ' "^ 100 for the symmetric form IT 



/3c^ 



^AT' 



which resembles the pattern 



74 



L% at N 



H 



Table g shows the first breaking points for the cold start P^ and for the hot start /3, 
These data are plotted in figure |3| for (3^ and figure ^ for (3^ . Again, we find clear linear 
dependence on A^ for /3^, as we found for the symmetric form. Additionally, we also find 
clear dependence on 1/A^ for /3^ . The fitted results are 
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0.0034A^ 
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Figure 2: Expectation value of tlie plaquette (top) and tlie Polyakov line (besides the top) versus 
the lattice coupling f3 ior N — 100 with the minimal skew-diagonal twist (cold start). As /3 is 



decreased, the ZiL symmetry is broken and restored as Z 
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64 
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0.2225 ±0.0025 


0.4625 ± 0.0025 


16 
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81 
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0.2125 ±0.0025 


0.5175 ±0.0025 
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0.2925 ± 0.0025 


0.3625 ± 0.0025 
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0.2075 ± 0.0025 


0.5875 ±0.0025 


36 
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0.2575 ± 0.0025 


0.3925 ± 0.0025 


121 


11 


0.2025 ±0.0025 


0.6525 ±0.0025 


49 


7 


0.2375 ± 0.0025 


0.4225 ± 0.0025 


144 


12 


0.1975 ±0.0025 


0.7325 ±0.0025 



Table 2: Critical lattice coupling from the weak coupling side /3^ and from strong coupling side 
13^ for the minimal skew-diagonal twist {m =1). 



where we used only A^ > 64 data for (5^, whereas all data are used for (3^ . As N is 

increased the /3^ approaches a point 0.190, where the phase transition Z^ — ^ Z^ takes 
place in the original EK model. These results suggest that the quantum fluctuation is so 
large that the Z^ symmetry is broken in exactly the same region as that in the original EK 
model. The lines for transitions (3^ and /3^ seem to intersect around the bulk transition 
point (3f ~ 0.35, which corresponds to A^ ~ 20 for the twist considered here. For smaller 
values than A^ ~ 20, we did not observe a signal of breakdown of the Z^ symmetry. 

Generic skew-diagonal twist 

Here, we show the numerical result of the generic twist ( p. 12 ). For this twist we use the 
skew-diagonal form because the Z^ symmetry breaking occurs at smaller A'^ than that in 
the symmetric form, which makes our investigation much easier. 
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Figure 3: Plot of (3^ versus N for the minimal 
skew-diagonal twist. The fit line is equation 



(3.5), which is obtained using iV > 64 data. 
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Figure 4: Plot of f3^ versus 1/A^ for the mini- 
mal skew-diagonal twist. The fit line is equation 



(3.6), which is obtained using N > 25 data. Ex- 
trapolation to 1/N = gives (]^ — > 0.18, which 
is close to the critical point in the original EK 
model, f3 = 0.19. 



We measure /3^ for this twist up to ?7i = 4. Table ^ shows the P^ for m = 2, 3, 4 and 
that for ?n, = 1 is presented in table |2[ These data are plotted in figure ^. From this figure 
we can find that the /3^ for each L are reduced as we increase m, and the dependence 
is linear in 1/m. The data at 1/m = in this plot are linearly extrapolated values. An 
interesting point is the behavior for the case L = 5. While the Z^ symmetry breaking 
is observed for ?7i = 1, 2 and 3, it is not seen for ?7i = 4 because the f3jf reaches the bulk 
transition point /3^ ~ 0.35 by increasing m. Figure ^ represents the same data in figure |5|, 
but the horizon axis is L^. As we have seen in the m = 1 case, the data for L > 8 are well 
fitted by the linear function of L^ for each m. From these figures, we find that the data 
for L > 8 are well fitted globally by a function: 



f3^^0mUL' + ^^ + 0.19. 
m 



(3.7) 



4. Discussions 

In this section we discuss the numerical results obtained in the previous section and the 
validity of taking the large- A^ and continuum limit for this model. 

4.1 Theoretical estimation of the Z^ symmetry breaking point 

In the previous section we showed our numerical results. In particular, we elaborately 
investigated /3^, the first Z^ breaking point from the cold start. From our investigation, 
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0.3525 ±0.0025 
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0.3675 ± 0.0025 
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0.3575 ±0.0025 
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0.3525 ± 0.0025 
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0.3925 ± 0.0025 
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0.3875 ± 0.0025 


196 


0.3825 ± 0.0025 
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0.4375 ± 0.0025 
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0.4275 ± 0.0025 


256 


0.4225 ± 0.0025 


9 
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0.4925 ± 0.0025 
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0.4825 ± 0.0025 


324 


0.4775 ± 0.0025 


10 


200 


0.5575 ±0.0025 


300 


0.5475 ± 0.0025 


400 


0.5425 ± 0.0025 



Table 3: /3^ for the generic skew-diagonal twist {m — 2,3,4). See also Tab. Hfor m = 1. 
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Figure 5: ^^ versus 1/m for L = 5, • • • ,10. Figure 6: Scaling of fi^ for m = 1,2,3 and 
/3^ for TO = oo is evaluated by extrapolating 4^ ^g ^Iso include /3^ for to = 00, which is 
these data with straight line. obtained by an extrapolation shown in figure ^ 



we found the clear linear behavior like ( |3.3| ), (3.5) and ( |3.7] ). These behaviors can be 
obtained through the following consideration. 



Energy difference between twist-eater r„ and identity Iat configurations 



We simply assume that the Z^ breaking is a transition from twist-eater phase C/^ = F^ 



Iat. For plainness, we consider Z"^ 



N 



% type 



to identity configuration phase C/^ 

breaking here. Of course we can treat Z^ — ^ Z^ — ^ Z^ — ^ Z]y — ^ Z^ (cascade) type 
breaking at a /3^, but the obtained behavior is not different from the former type. Firstly, 
we focus on the classical energy difference between these configurations. The energy can 



11 



be easily calculated from the action (|2.2|) as 

A5 = StekP^ = In] - StekP^, = r^] 

= /?iV25:|l-cos(^)U2vr^/3j;n^,. (4.1) 

For the generic twist, it becomes 

{247r^/3m^L^ (symmetric form), 

Stt"^ Pm'^ L'^ (skew-diagonal form). 

Note that the symmetric form is roughly three times more stable than the skew-diagonal 
form if both twists have equal quantum fluctuations. This is the reason that the Z^ 
symmetry breaking for the skew-diagonal form can occur at quite smaller A^ than that for 
the symmetric form, as is observed in our simulation. 

Quantum fluctuations and symmetry breaking 

Going away from the weak coupling limit, the system has quantum fluctuations. Here 
we naively expect that the Z^ symmetry is broken if the fluctuation around twist-eater 
configuration exceeds the energy difference AS. Because the system describes 0{N'^) 
interacting gluons, it is natural to assume that their quantum fluctuations provide 0{N'^) 
value to the effective action. For the generic twist the quantum fluctuation is 0{m?L'^). 
Combined with the fact ( |4.2| ), we can estimate the critical point f3^ as 

(5^ ~ L\ (4.3) 



which is consistent with the numerical results ( |3.3D , ( |3.5D and ([3.7|). In addition we can 
explain the difference of the coefficient of L^ in ( |3.3| ), ( |3.5D and ( |3.7| ) between the symmetric 
and the skew-diagonal form, which is roughly three times different, by the factor in (4.2). 



Although the above crude estimation reproduces the linear L behavior of /J^ ; we 
cannot explain the dependence on m. To catch the behavior completely, we need to make 
the discussion more sophisticated. However, we do not pursue this issue here because the 
m dependence can be negligible at the larger A^. 

This argument can be applied for other twist prescriptions like taking the twist phase 
as exp(z7r(L -I- \)/L), which is usually used for describing noncommutative spaces. (See 
appendix ^.) 

4.2 Continuum and large- A^ limit 

We have shown that the linear L^ dependence of the critical point (5^ could be explained 
by the theoretical discussion in this section. While our simulation is restricted in the small 
A^ region, we confirm that the behavior must continue to A^ = oo by combining with the 
discussion. Then the EK equivalence is valid only in the region /3 > /3^ ~ A^ even in the 



weak coupling limit and the large- A^ limit. As we mentioned in the section 2.4, both the 



ordinary YM with fixed physical volume and the NCYM theory with fixed noncommutative 
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parameter have essentially logarithm scaling ( 2.15 ) near the weak coupling limit. Then, 
because P^ grows faster than the logarithm, the EK equivalence does not hold in the 
continuum limit. 

5. Conclusions 

In order to study the nonperturbative nature of the large- A^ gauge theory by lattice simu- 
lations, the large- A^ reduction is very useful property for saving the computational effort. 
In this paper, we studied the phase structure of the TEK model, which has been a major 
way to realize the large-A^ reduction. Contrary to the naive hope in old days, at least 
in ordinary twist prescriptions as investigated in this paper, the Z^ symmetry is broken 
even in the weak coupling region and hence a continuum limit as the planar gauge theory 
cannot be described by the TEK model. For the NCYM, the situation is the same. We can 
also consider a lot of variation for the twist prescription and the combination of reduced 



and non-reduced dimension. For example, in |18, 19|, four-dimensional model with two 
commutative and two noncommutative directions was studied using two-dimensional lat- 
tice action. However, the Zjl^r symmetry is broken also in this model, and hence we cannot 
take a naive continuum limit. 

Another way for the reduction is the QEK model, in which the eigenvalues of the 
link variables are quenched. The QEK model might have no problem in principle, but its 
computational cost is larger than that of the TEK model. Although the TEK and QEK 
model are reduced models to a single hyper-cube, recent studies deviate from them. The 



contemporary method might be the partial reduction |22|. This work showed that the 
large-N reduction is valid above some critical physical size Ic- This means that for a lattice 
size L the reduction holds below some lattice coupling P{L). In order to take continuum 
limits we should avoid the bulk transition point P^ , causing the condition /3^ < (3{L) to 
be necessary. That is, there is a lower limitation for the lattice size L^ for the continuum 
reduction. In addition, the twist prescription is also applicable to the partial reduction pSJ . 
Due to the twisted boundary condition, the lower limitation Lc can be reduced. Therefore, 
combination of the twist prescription and the partial reduction would be quite efficient in 
the current situation. 

Note also that NCYM on fuzzy R^ could be realized by using TEK with quotient 
conditions [^ which give a periodic condition to eigenvalues and hence quantum fluctuation 
is suppressed. Further study in this direction would be important. 

Acknowledgments 

The numerical computations in this work were in part carried out at the Yukawa In- 
stitute Computer Facility. The authors would like to thank Sinya Aoki, Hikaru Kawai, 
Jun Nishimura, Masanori Okawa, Yoshiaki Susaki, Hiroshi Suzuki and Shinichiro Yamato 
for stimulating discussions and comments. M. H. was supported by Special Postdoctoral 
Researchers Program at RIKEN. T. H. would like to thank the Japan Society for the 
Promotion of Science for financial support. 



13 



A. Double scaling limit as the noncommutative Yang-Mills theory 

The TEK model can be used to formulate gauge theories on noncommutative spaces non- 



perturbatively [0, ^, 24|. In this appendix, we give a review for the construction of the 
NCYM from the TEK model Q, a discussion for the scaling and some supplemental com- 
ments for our analysis. 

By taking f/^ = e* f*, where a corresponds to the lattice spacing, and expanding the 
action of the TEK model ( |2.2| ), we have its continuum version 



STEK,continuum = ~2~2 \^ '^^ ([^M; ^A ~ ^(^ fiv) (^.1) 



up to higher order terms in a, where 

Then, by expanding the action around a classical solution of ( [A.l| ) 

4°^=Pm' [P^.Pu\=iO^,u, (A.3) 

we obtain the U{1) NCYM on fuzzy M"^ as follows. Let us define the "noncommutative 
coordinate" x^ = (^~^) Pu- Then we have 

\^j;i^^x''] = -i{d-^Yr (A.4) 

This commutation relation is the same as that of coordinate on fuzzy R^ with noncommu- 
tativity parameter 9, and hence functions of x can be mapped to functions on fuzzy M^. 
More precisely, we have the following mapping rule: 



(A.5) 



fix) = Ek /(^)e*'^ -- /(^) = Ek /(fc)e*'", 
f{x)g{x) ^ f{x)'kg{x), 

where * represents the noncommutative star product, 

f{x)*g{x) = f{x)exp(^-'-l)^{e-'r''l),^ g{x), (A.6) 

and we obtain C/(l) NCYM action 



Suii)NC = -7^— d x F^^-kF^^ (A.7) 

with coupling constant 

glc = 47r29V^/dit^. (A.8) 

In the same way, by expanding the action (|A.1|) around Ajj^ = p^® Im, U{m) NCYM can 
be obtained. From ( p. 13 ), it is apparent that the generic twist gives the U{m) NCYM. 
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Intuitively, the vacuum configuration ( |2.13D describes m-coincident fuzzy tori and fuzzy M^ 
is realized as a tangent space. 

In order to keep the nonconimutative scale finite, we should take the double scaling 
limit with 

One-loop beta function for U{m) NCYM is given by 



2 11 o ^ 

/?i-ioop (s-ivc) = ~7^^^9nc + 0{gNc)- (A.IO) 

Therefore, the 't Hooft coupling (3 scales as 

/3^^ logiV. (A.ll) 

9nc 

Then, the scaling we should take for the NCYM is just the same as that for the ordinary 
YM, and Z^ symmetry is broken in the scaling limit. Therefore, fuzzy torus crunches to 
a point and hence the fuzzy M.^ cannot be realized^. 

Of course, we can also use other twist prescriptions. In order to make the periodicity 
of the discretized fuzzy torus correct, we usually take the twist as exp(i7r(L + 1)/L) Q. 
Regardless of the difference of the twist, the conclusion might not be altered. Here we 



repeat the discussion in section |4.l| . In this case, the Z^ is likely to break down to Zg. The 
difference between potentials in twist-eater and Z^-preserving configurations is 



A5 ~/3iV2{l- cos ('yU ~/3m2L2, (A.12) 



which is the same order as ( [4.1| ). Then the behavior of the critical point /3^ ( |4.3| ) is not 
changed. 
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